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Deflection angle for a light ray travelling in the equatorial plane of a rotating Kerr mass has been 
calculated in past by various investigators. Considering the light ray to be travelling only slightly 
above the equatorial plane, calculations have been made in the present paper for such a ray for its 
deflection angle. We calculate deflection angles for the light ray at various heights, which are small 
compared to the impact parameter and derive corresponding analytical expressions for deflection 
angle. 


I. INTRODUCTION 

In general relativity, Schwarzschild line element gives 
the solution of Einstein’s field equation for the exterior 
of an uncharged, static mass [1]. Bending angle for 
Schwarzschild field was calculated by many authors up 
to second and higher order terms. Keeton et al. [2] cal¬ 
culated the light deflection angle for Schwarzschild black 
hole in weak field approximation. Iyer and Petters [3] cal¬ 
culated light deflection angle for strong field conditions, 
which under weak field approximation is in exact match 
with that of Keeton et al. [2]. 

The solution of Einstein’s field equation which describes 
the exterior of an uncharged, rotating mass was obtained 
by Kerr in the year of 1965 [4], According to this line 
element, the frame dragging an unusual prediction of the 
general relativity is exhibited by such rotating masses. 
The prediction of this effect is that, all objects coming 
close to a rotating mass would be entrained to participate 
in its rotation, not because of any applied force or torque 
that can be felt, but because of the curvature of space 
time associated with the rotating mass. At close enough 
distance all objects even light must rotate with the mass. 
So the light deflection angles for Schwarzschild and the 
Kerr geometries are not the same. There are some extra 
terms, which arise due to the effect of rotation. 

Bending angle for Kerr field in equatorial plane (i.e 
9 = was calculated by Iyer et al. [5, 6]. Accord¬ 
ing to their result, deflection produced in presence of a 
rotating black hole explicitly depends on direction of mo¬ 
tion of the light. If the light ray is moving in the direc¬ 
tion of rotation (called the proggrade orbit), then the 
deflection angle is higher than that for the zero rotation 
Schwarzschild field. If the light ray is moving opposite to 
the direction of rotation (called the retrograde orbit), the 
bending angle will be smaller than that for Schwarzschild 
field. Dubey and Sen [7] had used Kerr geometry, to show 
how the gravitational redshifts are affected as photon is 
emitted at various latitudes of rotating mass. Kraniotis 
[8] derived the analytic solutions for the lens equations in 
terms of Appell and Lauricella hypergeometric functions 
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and the Weierstrafi modular form for a Kerr mass. All 
the above mentioned calculations (for Schwarzschild and 
Kerr) were done using the null geodesic for light ray. 

On the other hand, another technique has been devel¬ 
oping by some authors called effective refractive index of 
material medium. Here the effect of gravitation is seen as 
the change in the refractive index of the medium through 
which light is travelling. In the year 1958, Balaz [9] used 
this method to calculate the effect of gravity due to a 
rotating mass on the direction of polarization vector of 
an electromagentic wave. Atkinson [10] used this method 
to study the trajectory of light ray near a very massive, 
static and spherically symmetric star [10]. Fischbacli and 
Freeman [11] derived the second order term of gravita¬ 
tional deflection using the same method. Sen [12] and 
Roy and Sen [13] calculated the light deflection angles 
for Schwarzschild and Kerr masses using this method. 
All these above mentioned works for Kerr field were 
done considering the light ray contained in the equatorial 
plane. But there are authors who have worked on non- 
equatorial plane also. Bozza [14] obtained the lensing 
formula, and calculated the relativistic image position 
for a light ray trajectory close to equatorial plane of a 
Kerr black hole. However, he did not calculate the light 
deflection angle explicitly for the non-equatorial plane. 
Aazami et al. [15, 16] worked on quasi-equatorial regime 
of Kerr black hole and calculated the two components 
of light bending angle, along the direction of equatorial 
plane and perpendicular to the equatorial plane of a Kerr 
black hole. Their expression for light deflection angle 
along the direction of equatorial plane is in exact match 
with the result of Iyer et al. [6]. 

For a static charged gravitational field, the solution of 
Einstein field equation was first given by Reissner [17] 
and Nordstrom [18] independently, known as Reissner- 
Nordstrom solution. Eiroa et al [19] calculated the light 
deflection angle for this line element and showed that 
charge itself has some effect on curvature of space-time. 
Vibhadra et al. [20] obtained the light bending angle 
for Janis-Newman-Winicour space-time (another solu¬ 
tion for static, charge mass) and it was shown that for 
zero static charge, it reduces to Schwarzschild bending 
angle. The exact solution of Einstein’s field equation for 
a rotating charged mass was found by Newman et. al 
[21, 22] in the year 1967, which is known as Kerr-Newman 
metric. Chakraborty and Sen [23] have recently obtained 



2 


Spin Axis 



p 



Equatorial Plane 


FIG. 1. Geometry of the equatorial plane of a rotating mass, 
where the light ray is passing through the point P, perpendic¬ 
ular to the plane of paper. Here, l is the height of the light 
source above the equatorial plane, u is the projected impact 
parameter on the equatorial plane and r o is the closest ap¬ 
proach of the light ray. 

the light deflection angle for a charged, rotating mass in tational body, where gravity dominates over the rotation, 

the equatorial plane. Under the condition of zero charge, So the light ray remains on the gravitational plane i.e. all 

their expression for deflection angle reduces exactly to over its trajectory light ray maintains a height l above the 

that of Iyer et al. [6]. equatorial plane. Direct consequence of this assumption 

In this paper we consider an uncharged, rotating Kerr is that the latitude angle 9 will be a constant. Strictly 

mass where the light ray is coming from an off-equatorial speaking as the light ray undergoes deflection the height 

source i.e. the source at infinity is not situated on the will not maintained constant all over its trajectory. How- 

equatorial plane, but it is placed at a very small height l ever, for the ease of calculation we assumed the height to 

above the equatorial plane (Fig 1). For a rotating body, be constant and such approached have been also adopted 

rotation drags the trajectory of light ray towards the ro- by other authors in past. With this assumption we calcu- 

tation plane i.e. the equatorial plane where the rotation late in the present work the off-equatorial light deflection 

has its maximum effect. angle for a slowly rotating body. 

In this paper we consider a slowly rotating, massive gravi- 


II. NULL GEODESIC IN KERR FIELD 


where 


The Kerr line element in Boyer-Lindquist co-ordinate 
(r, 8, (f), t) can be written as [24, page 313] 


2 i 2 2/i 

p = r + a cos 0 


( 2 ) 


A = a 2 +r 2 — 2 mr 


( 3 ) 



and to = and a = further c is the velocity 
of light in free space, M is the mass of the gravitating 
body, G is the gravitational constant and J is the angular 
momentum of the gravitating body. Lagrangian (£) can 


2 mra 2 sin 2 8 



cdtdcj) (1) 
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be written as, [24, page 96] 

£ = (4) 

where the dot indicates derivative with respect to some 
affine parameter. So the lagrangian (C) can be written 
as [24], 


„„ 2mr, 9 -o Aamr sin 2 9 . • p 2 • 9 •„ , 9 , 

2 £ = (1 -—)c 2 f 2 +--- ct(p— h —r 2 —p 2 9 2 —(r 2 +a 2 

p 1 p z A 


+ 



) sin 2 9(f> 2 


(5) 


From this equation, the energy (E) and angular momen¬ 
tum in the direction of rotation (L,) can be written as 
[24, page 346]: 


E = P,= 


dC 

d(ct) 


2 Mr. • 
= (1 - 2~)ct 


2 amr sin 6 • 


( 6 ) 


dC 2amr sin 2 0 . 

= —P<t> = JTT = 5-cf+ 

d(p p 2 


2 a 2 mr sin 2 


9 


sin 2 9(f) 


(7) 


It is possible to show that impact parameter £ s = s£ = 
[24, page 123]. Here, s defines the direction of mo¬ 
tion of the light ray. For s = +, the light ray is in the 
direction of rotation of the mass and for s = —, the light 
ray is in the opposite direction of rotation [ 6 ]. Here 
is the first constant of motion. The second constant of 
motion is p = where \ = K — ( L z — aE ) 2 . The quan¬ 
tity x is known as the Carter constant [25] and K is the 
conserved quantity along the direction 9. 

Following Carter [24, 25], we can have the four geodesic 
equations governing the motion of light ray in Kerr ge¬ 
ometry as: 


[2mra + (p 2 — 2rnr)£ esc 2 9]E 
(a 2 + r 2 — 2 mr)p 2 


( 8 ) 


[(r 2 + a 2 ) 2 — Aa 2 sin 2 9 — 2 mra£,]E 

ct =-A7- (9) 


(77 + a 2 cos 0 — £ 2 cot 9) 2 E 


( 10 ) 


[r 4 — (77 + £ 2 — a 2 )r 2 + 2mr(a — £ s ) 2 + 2mrr] — a 2 ij] 2 E 
r P 2 

( 11 ) 

Now let us consider that the source of the light ray is not 
contained in the equatorial plane, but slightly above it 


by a distance l. Under this condition 9 is not equal to -|. 
If there is no gravitational held, the light ray will follow 
a straight line. The projection of this straight line on 
equatorial plane will be at a distance u from the origin of 
the rotating mass (Fig 1), which we can call the projected 
impact parameter. At this minimum projected distance, 
the light ray has a height l from the equatorial plane and 
?/> is the angle that the light ray makes with the equatorial 
plane [14]. So, 

tan ip = — ( 12 ) 

u 

The value of ip is very small, as l is only slightly above 
the equatorial plane. In other words l is much smaller 
then u (l C«). 

These two parameters (l,u) can be used to level the light 
ray coming from infinity. The relation between the two 
parameters (l,u) and two constants of motion (£ s and rf) 
can be written as [14] : 

= u cos ip (13) 

and 

77 = l 2 cos 2 ip + (u 2 — a 2 ) sin 2 ip (14) 


Now the relation between colatitude angle 9 and the angle 
ip can be written as, 

9=^-ip (15) 


On the equatorial plane (i.e. I = ip = 0), = u and 

77 = 0 . 

Now, the geodesic equations can be written in terms of 
ip by using equation (15) in equation ( 8 , 9, 10, 11) as 
follows: 


Or, 


■ [2 mr + (r 2 + a 2 sin 2 ip — 2mr)^ s sec 2 ip\E 

( a 2 + r 2 — 2mr)(r 2 + a 2 sin 2 ip) 

[2 mra + r(r — 2m)^ s sec 2 ip + a 2 £ s tan 2 ip\E 
(a 2 + r 2 — 2mr)(r 2 + a 2 sin 2 ip) 


(16) 


7 /) 


(77 + a 2 sin ip — £ 2 tan ip) 2 E 
r 2 + a 2 sin 2 ip 


(17) 


[(r 2 + a 2 ) 2 — Aa 2 cos 2 ip — 2mra^ s \E 
A(r 2 + a 2 sin 2 ip) 

[ 7 ’ 4 — (77 + £ 2 — a 2 )r 2 + 2mr(a — £ s ) 2 + 2mrr) — a 2 i 7 ] 5 E 
r 2 + a 2 sin 2 ip 

( 19) 

Here, we considered ip is very small i.e. ICti. So we can 
expand the trigonometric functions in the above geodesic 
equations up to second order terms as: 

tan 2 ip ~ sin 2 ip ~ ip 2 


( 20 ) 
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sec 2 ip ~ 1 + ip (21) 

we can have the new format of geodesic equations by 
using equations (20), (21) in equations (16) to (19) as: 


It may be noted here that, ro = y/u 2 + l 2 ■ Further from 
equation (25), we can write : 

r o - (V + £« - a 2 )rl + 2 mr 0 (a - £ s ) 2 + 2 mr 0 ij - a 2 rj = 0 


[2mra + r(r — 2m)£ s (l + ip 2 ) + a 2 £ s ip 2 ]E 
( a 2 + r 2 — 2 mr)(r 2 + a 2 ip 2 ) 


a 2 2 to a 2 r] 


= (!-«)-— ( 1 -^) 


6 ' £ 


(?7 + a 2 ip — £ 2 ip) 2 E 


2 mij rja 2 
Ms + ts r o 


• [(r 2 + a 2 ) 2 - Aa 2 (l - V> 2 ) - 2mra^ s ]£’ 

C '~ A(r 2 + a 2 ip 2 ) ( j 


III. DEFLECTION ANGLE FOR OFF 
EQUATORIAL SOURCE 


Under the above kind of geometry, one may write the 
. _ [r 4 — (77 + £ 2 — a 2 )r 2 + 2mr(a — £ s ) 2 + 2mrr] — a 2 r)\ 2 E expression for the light deflection angle as : [26, page 

r 2 + a 2 ip 2 189] 

(25) 

In the present case, one can show that, r obtains a local = 9 / ^ (27) 

extremum for the closest approach r 0 , so that we can J ro dr 

write: 

Using equation (22) and (25) in equation (27), one may 
write 


r\ r =ro — 0 


A<p = 2 


[2 mra + r(r — 2 m)£ s (l + ip 2 ) + a 2 t; a ip 2 ]dr 
[a 2 + r(r — 2 m)][r 4 — (77 + £ 2 — a 2 )r 2 + 2 mr(a — £ s ) 2 + 2mi'ij — a 2 rj]i 


A(P = 2 


r°° M=+ r (r-2m.){l + ip 2 ) + a 2 ip 2 ]dr 

r ° [a 2 + r(r - 2m)]£ s r 2 [£ - ^ + ^(1 - £) 2 - + y%(2 mr - a 2 )]^ 


As was done by some previous authors [15, 16], here we substitute G = 1 — (^-) 2 = 1 — a 2 (|j) 2 and F = 1 — (-^) = 
1 — sa where a = A. So the new form of the above deflection expression becomes: 


A 0 = 2 


_ [1- gfU + V> 2 (l- gf + ^)]dr _ 

[a 2 + r(r - 2 m)] [ 4 - - # + + 73 W ( 2 mr - a 2 )]^ 


We substitute, h = x = ^ and n= Accordingly we have, 


dx = —^dr 
r 2 


and the limits will now change, when r —> 00 , we have x — > 0 and when r —> ro, then x — > 1. Using these 
substitutions, in above equation we have: 


A</> = 2 


f 1 [1 — 2hxF + ip 2 (l — 2hx + a 2 h 2 x 2 )\dx 

0 [1 — 2/ix + d 2 /i 2 ;r 2 ][4§ — Gx 2 + 2hF 2 x 3 — nx 2 + 2 hnx 3 — na 2 h 2 x 4 ]? 


Using the above substitutions, we can rewrite the equation (26) as : 


rn a 2 . 2 m., a .0 ?7 2m?7 77a 2 

t | = 1 - 7 ^ — 1 +72-—d + 

ts Zs r 0 6 Cl U>£f £ s 2 ?o 
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= G — 2 hF 2 + n — 2/m + na 2 h 2 
Combining equation (31) and (32), we can write 

[1 — 2 hxF + '0 2 (1 — 2 hx + d 2 h 2 x 2 )]dx 


A(j) = 2 


or, 


o [1 — 2/ix + d 2 /i 2 x 2 ][G(l — a: 2 ) — 2hF 2 (l — x 3 ) + n(l — x 2 ) — 2hn(l — x 3 ) + na 2 h 2 ( 1 — a: 4 )] 2 

[1 — 2hxF + ^> 2 (1 — 2/ix + a 2 h 2 x 2 )]dx 


Acj) = 2 


(32) 


(33) 


(34) 


/0 [1 — 2hx + a 2 h 2 x 2 ][(G + n)(l — a; 2 ) — 2 h(F 2 + n)(l — a; 3 ) + na 2 h 2 ( 1 — a: 4 )] 2 

Let us assume G + n = g and F 2 + n = f. If the light ray is contained in the equatorial plane, then we shall have 
n — 0 and we have, G = g and f = F 2 . So the new form of the equation (34) will be, 

[1 — 2hxF + ip 2 ( 1 — 2 hx + a 2 h 2 x 2 )\dx 


or, 


or, 


or, 


A(j> = 2 


A(j> = 2 


^ Jo [1 — 2hx + a 2 h 2 x 2 ][g(l — x 2 ) — 2/i/(l — a; 3 ) + na 2 h 2 (l — a; 4 )]2 

^ j" 1 [1 — 2hxF + ip 2 {\ — 2hx + d 2 h 2 x 2 )]dx 

Jo [1 - 2 hx + d 2 h 2 x 2 ]\Jg{\ - x 2 )[l - + x 2 )]? 

y 1 [1 — 2hxF + ip 2 {\ — 2hx + a 2 h 2 x 2 )\dx 

'° [1 - 2 hx + d 2 h 2 x 2 } ^g(l-x 2 )^l~ 2hj^[l + 2^(1 + x 2 ){l - 2 

[ — - [1 — 2hxF + f/; 2 (l — 2 hx + d 2 h 2 x 2 )][\ — 2 hx + d 2 h 2 x 2 F 1 { 1 — 2/j^-2 — 

f 0 \ J 5(1 — x 2 ) 5(1 - X 2 ) 

"" 2 h 2 9,.. /(1-a; 3 ) ! 1 


(35) 


[1 H-(1 + x ){1 — 2 h- 


-r 1 ] 5 - 


9 ~ 'g( I-* 2 ) 

For weak deflection limit, following [15, 16] one can assume, m <C r 0 , in other words, /i <C 1. So the above equation 
can be expanded in Taylor series in terms of h. By considering terms up to second order only, we can write : 


A<p = 2 


dx 


0 x/fi'C 1 - x 2 ) 


[1 — 2 hxF + f/r(l — 2 hx + ah 2 ar)][l + 2hx + ar h “(4 — a 2 )} 


fh 1 — a; 3 3/ 2 /i 2 1 — a; 3 nd 2 K 2 


[! + —( 1 -^ 2 ) 

q 1 — x z 


-(- 




g L — x*' 2 g 2 y l-x 2 ' il 2g 

Multiplying term by term and retaining only up to second order of h we get, 


(1 + x 2 )] - 7 r 


(36) 


A<j> = 2 


dx 


0 \fg( 1 - x 2 ) 


f (1 _ o-S'l 

[1 + /i{2x(l — F) + ^ 


2 n 2 _ 1 2/x(l — F) (1 — x 3 ) 3/ 2 (1 — x 3 ) 2 2 \ na 2 2 w < 2 mi 

(4 - “ - 4F > + —- + v (Tv^* 1 + - 1 Y 1 + 1 )(1 + * -' 


(37) 


Now integrating term by term and retaining terms only up to second order in h 1 we can write : 

A* = - 1) + 4 M /(1 + ’ / ’ 2 > + 9 - Fg | + )>*[-4{/(l + 1?) + 9 - 


V5 ff 2 52 

“7~~~-t{ 15/ 2 (1 + V> 2 ) - 4ff(F - 1)(3/ + 2g) - 2g 2 d 2 } - (1 + ^ 2 )] 

4 15^2 4^2 


(38) 


The above equation (38) represents the off-equatorial de- This expression is a function of mass, rotation parameter 
flection of light due to a rotating mass in weak field limit. 
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and ip only. We verify our result under some limiting 
conditions below: 

If the light ray is contained in equatorial plane only, then 
ip, n will be equal to zero and g = G, f = F 2 . Under 
these conditions we find, the above equation (38) will 
reduce to the following, 


= 7r (y=- 1 )+ 4 ^[ 


F+G-FG 
Gi 


]+h 2 [—4 


F{F+G-GF) 

Gi 


+ -r-77^{^F 2 -4G(F-l)(3F+2G)-2G 2 a 2 }\ (39) 
4 15G 2 

The above equation (39) is same as the expression for 
equatorial deflection of light by Kerr mass obtained by 
Aazarni et al. [16, equation no. (B17)]. 


IV. DISCUSSION AND CONCLUSION 

The main focus of the paper is to study the light de¬ 
flection angle (A <p) as a function of the height ( l ) of the 
light source above the equatorial plane in Kerr geome¬ 
try. As obtained in equation (38), the deflection angle 
(A <p) is a function of ip ~ where u is the projected 
impact parameter. So by varying the source height ( l ) 
and projected impact parameter (u), the main purpose 
of the paper can be achieved, to study their effect on the 
deflection angle. 

In order to study the relation between A tp and l we obtain 
the values of Atp as a function of tan ip = ^ in Tablel, 
by considering, the radius of Sun 6.955 x 10 8 meter as 
the closest approach of light ray (ro) for both prograde 
and retrograde motions. The same was repeated, by con¬ 
sidering a slow rotating pulsar PSR J 1748-2446 [27] for 
prograde and retrograde orbits and the results are repro¬ 
duced in Fig.2. For these calculations, for the pulsar the 
equatorial projected impact parameter ( u ) is taken to be 
u = 7500 x 2m, 7501 x 2m, 7502 x 2?n to maintain the 
weak field approximation. Subsequently, to understand 
the relation between A (p and u, the plots were repeated 
in all the above mentioned cases, for three different val¬ 
ues of u/2m. Here we should mention that, the pulsar 
has been chosen completely arbitrarily, except that the 


ratio of Kerr parameter(a) to Schwarzschild radius(2m) 
is same as that of Sun approximately, which will help us 
to compare results later. For that reason we choose a 
realistic pulsar PSR J 1748-2446. 

In Fig2, we draw the deflection values (A<^>) against iemip 
for the two cases (prograde and retrograde) for different 
parametric values of u/2m.. 

By studying the nature of the graphs, it is clear that the 
light deflection angle increases with the r/>(~ -) for both 
prograge and retrograde motion. For Sun (Tablel) and 
pulsar (Fig2) the above variation of A <p as a function of 
ip is clearly seen. On the other hand, it is quite obvious 
that the deflection angle decreases with the increase of 
projected impact parameter normalized by Schwarzschild 
radius [u/2m) (cf. Fig2 and Tablel). 

All the above mentioned calculations were done by con¬ 
sidering ip to be very small (l <C u). Since from our 
expression we know that A (p is a function of both r and 
ip, therefore if we consider l to be in the same order as u, 
then instead of ip being a constant, it would be a func¬ 
tion of r itself and then we have to write the expression 
of light deflection angle as: 


or, 


A<p = 2 


dg>{r,ip{r)} 

dr 


.dr — 7r 


A<P = 



dtp dip 
dip' dr 


— 7 T 


We may try to deduce the the values of and 
utilizing equation (16), (17) and (19). However, this we 
postpone for our future work and in the present work we 
report our work under the condition of l <C u only. 
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TABLE I. Variation of bending angle of light (A <j>) as a function of ■)/> ~ 
Sun. To note, at the location of solar radius, value of u/2m = 2.35 x 10 s 

— for three different values of u/2m with a = 0.5 for 
at ip = 0. 

u/2m 

1>~± 

T U 

bending angle for prograde motion (arcsec) bending angle for retrograde motion (arcsec) 

2.35 x 10° 

1 x 10“° 

1.753929 

1.753933 


2 x 10 -5 

1.754125 

1.754129 


3 x 10 -5 

1.754451 

1.754455 


4 x 10 -5 

1.754906 

1.754910 

2.82 x KV 

1 x 10“° 

1.461618 

1.461621 


2 x 10“ 5 

1.461814 

1.461817 


3 x 10“ 5 

1.462140 

1.462142 


4 x 10“ 5 

1.462595 

1.462597 

3.29 x 10° 

1 x 10“° 

1.252825 

1.252826 


2 x 10 -5 

1.253021 

1.253023 


3 x 10“ 5 

1.253346 

1.253803 


4 x 10 -5 

1.253801 

1.253803 


3.29 x 10' 






FIG. 2. Bending angle (arcsec) as a function of tan. ip = ^ 
with constant rotation (a = 0.5) for both prograde and ret¬ 
rograde motion of PSR J 1748-2446 for different values of 
projected impact parameter (u/2m). 
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